Fluid dynamics and noise in bacterial cell-cell and cell-surface scattering 
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Bacterial processes ranging from gene expression to motility and biofilm formation are constantly 
challenged by internal and external noise. While the importance of stochastic fluctuations has been 
appreciated for chemotax;is, it is currently believed that deterministic long-range fluid dynamical 
effects govern cell-cell and cell-surface scattering - the elementary events that lead to swarming 
and collective swimming in active suspensions and to the formation of biofllms. Here, we report the 
first direct measurements of the bacterial flow fleld generated by individual swimming Escherichia 
coli both far from and near to a solid surface. These experiments allowed us to examine the relative 
importance of fluid dynamics and rotational diffusion for bacteria. For cell-cell interactions it is 
shown that thermal and intrinsic stochasticity drown the effects of long-range fluid dynamics, im- 
plying that physical interactions between bacteria are determined by steric collisions and near-fleld 
lubrication forces. This dominance of short-range forces closely links collective motion in bacterial 
suspensions to self-organization in driven granular systems, assemblages of biofilaments, and animal 
flocks. For the scattering of bacteria with surfaces, long-range fluid dynamical interactions are also 
shown to be negligible before collisions; however, once the bacterium swims along the surface within 
a few microns after an aligning collision, hydrodynamic effects can contribute to the experimentally 
observed, long residence times. As these results are based on purely mechanical properties, they 
apply to a wide range of microorganisms. 
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Collective behavior of bacteria, such as biofilni for- 
mation [l|, swarming [2| and turbulence-like motion in 
concentrated suspensions [3|, |j] , has profound effects on 
foraging, signaling, and transport of metabolites [sl, |6[, 
and can be of great biomedical importance [7|, |8[ . Large- 
scale coherence in bacterial systems typically arises from 
a combination of biochemical signaling [9| and physical 
interactions. Recent theoretical models that focus on 
physical aspects of bacterial dynamics idei itify pairwise 
long-range hydrodynamic interactions J10l - ll6l | as a key 
ingredient for collective swimming. Such "microscopic" 
approaches underpin continuum theories that aim to de- 
scribe the phenomenology of microbial suspensions |17l - 
|23| . An assumption underlying many of these theories is 
that a self-propelled bacterium can be modeled as a force 
dipolc; its body exerts a drag force F on the fluid that 
is balanced by the rearward flagellar thrust —F. The 
leading order fluid velocity fleld at distance r is therefore 
a dipolar "pusher" flow of magnitude u oc Fi/rjr'^ (see 
streamlines in Fig. [TJB), where rj is the viscosity, and i 
the distance between the forces [2J, |25[ . While higher 
order corrections may be due to force-quadrupole contri- 
butions [26[, the hypothesis that the leading-order flow 
field around a bacterium is dipolar has not yet been ver- 
ified experimentally. 



A closely related, controversially discussed issue [27- 
|30| is the relevance of long-range hydrodynamic interac- 
tions in the scattering of bacteria with surfaces, a phe- 
nomenon intimately linked with surface accumulation 



and biofilm formation. Cell-surface scattering is very 
similar to cell-cell scattering, since, by analogy with im- 
age charges in electrostatics, a bacterium that swims near 
a surface induces an "image bacterium" on the opposite 
side of the wall to yield the no-slip boundary condition 
on the surface [27|, |31[ ; bacterium-surface scattering can 
therefore be analyzed as the interaction of a bacterium 
with its hydrodynamic image. Several recent calculations 
for microswimmers near surfaces have found that pusher- 
type organisms (those with thrust generated behind the 
cell body) should display a passive stable alignment of the 
swimming direction with the wall |27l. |32h35|. However, 
direct measurements of the three-dimensional swimming 
tracks of bacteria near surfaces suggest that they simply 
collide with the surface [28J, |36[ . This raises the question: 
Is this discrepancy between experiment and theory due 
to incomplete knowledge of the bacterial fiow field, or is 
the magnitude of the flow so small that fluid-mediated 
interactions are irrelevant? 

The need for experimental tests of the force dipole 
assumption and, more generally, of the relevance of 
fluid-mediated interactions for bacteria, is further illus- 
trated by recent measurements on Chlamydomonas rein- 
hardtii [37J, l38[, the archetypal "puller" microorganism 
(thrust generation in front of the cell body) that was 
thought to generate a simple force dipole flow with op- 
posite sign to the bacterial one [25[- Surprisingly, these 
experiments showed that while such a dipolar flow exists 
at large distances from the organism, in regions where 
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FIG. 1: Average flow field created by a single freely swimming bacterium far from surfaces (A-D), and close to a wall (E- 
H). Streamlines indicate the local direction of flow, and the logarithmic color scheme indicates flow speed magnitudes. (A) 
Experimentally measured flow fleld in the bacterial swimming plane, with the inset showing the anterior-posterior asymmetry 
close to the cell body. (B) Best-flt force dipole flow. (C) Residual flow fleld, obtained by subtracting the best-fit dipole model 
from the measured field. (D) Radial decay of the flow speed u in different directions, with r — corresponding to the center of 
the cell body. For distances r < 6 /im the dipole model overestimates the flow field behind and to the side of the cell body. (E) 
Experimentally measured fiow field in the bacterial swimming plane, for bacteria swimming parallel to a wall at a distance of 2 
/xm. (F) Best-fit force dipole fiow, where the presence of the wall causes inward and outward streamlines to join. (G) Residual 
flow field. (H) The fiow speed decays much faster for bacteria swimming close to a wall, as the fluid velocity must vanish on 
the surface. 



the flow magnitudes are significant (more than 1% of the 
swimming speed), the flow topology is qualitatively dif- 
ferent, and more accurately described in terms of a triplet 
of force singularities (one for the cell body and one for 
each flagellum) [37| . Here we report the first direct mea- 
surement of the flow fleld around individual freely swim- 
ming bacteria, using Escherichia coli as a model. We find 
that the pusher force dipole provides a good approxima- 
tion to the flow field both when the organism is far from 
surfaces and close to a no-slip boundary, yet the magni- 
tude of the flow is very low. Using the experimentally 
determined flow field parameters, the hydrodynamic in- 
teraction of two E. coli can be calculated, and it is found 
to be washed out by rotational diffusion of the swimming 
direction for closest encounter distances > 3 /im ~ a re- 
sult that should hold for many other bacterial species 
due to the similarity of motility parameters. Similarly, 
analysis of cell-surface encounters suggests that hydrody- 
namics plays a negligible role, except when a bacterium 
swims along a surface at a small distance (less than a few 
microns) after an inelastic aligning collision. In this case, 
hydrodynamic effects can contribute to the observed long 
residence times near surfaces. 



RESULTS 
Flow Field Far from Surfaces 

To resolve the minisculc flow field created by bacteria, 
individual gfp-\ahe\ed, non-tumbling E. coli were tracked 
as they swam through a suspension of fluorescent tracer 
particles (see Materials and Methods). Measurements far 
from walls were obtained by focusing on a plane 50 ^m 
from the top and bottom surfaces of the sample cham- 
ber, and recording '^ 2 terabytes of movie data. Within 
this data we identifled ~ 10^ rare events when cells swam 
within the depth of fleld (2 /zm thick) for > 1.5 s. By 
tracking the fluid tracers during each of the rare events, 
relating their position and velocity to the position and 
orientation of the bacterium, and performing an ensem- 
ble average over all bacteria, the time-averaged flow field 
in the swimming plane was determined down to 0.1% of 
the mean swimming speed Vq = 22 ± 5 ^m/s. As E. 
coli rotate about their swimming direction, their time- 
averaged flow field in three dimensions is cylindrically 
symmetric. The present measurements capture all com- 



ponents of this cylindrically symmetric flow except the 
azimuthal flow due to the rotation of the cefl about its 
body axis. In contrast with the flow around higher organ- 
isms such as Chlamydomonas |37l.l38l| and Volvox [37[, the 
topology of the measured bacterial flow field (Fig. lA) 
is that of a force dipole (shown in Fig. IB). Yet, there 
are some differences between the force dipole fiow and 
the measurements close to the cell body, as shown by the 
residual of the fit (Fig. IC). 

The decay of the fiow speed with distance r from the 
center of the cell body (Fig. ID) illustrates that the mea- 
sured flow field displays the characteristic 1/?'^ form of 
a force dipole. However, the force dipole model signifi- 
cantly overestimates the fiow to the side and behind the 
cell body, where the measured flow magnitude is nearly 
constant over the length of the flagellar bundle. The force 
dipole fit to the far field (r > 8 fim) was achieved with 
two opposite force monopoles (Stokeslets) at variable lo- 
cations along the swimming direction. As r = corre- 
sponds to the center of the cell body in Fig. ID, and not 
the half-way point between the two opposite Stokeslets, 
the fit captures some of the anterior-posterior asymmetry 
in the flow magnitude u. From the best fit, which is in- 
sensitive to the specific algorithms used, we obtained the 
dipole length £ = 1.9 yum and dipole force F = 0.42 pN. 
This value of F is consistent with optical trap measure- 
ments [39| and resistive force theory calculations [40| . It 
is interesting to note that in the best fit, the cell drag 
Stokeslet is located 0.1 /zm behind the center of the cell 
body, possibly refiecting the fluid drag on the flagellar 
bundle. 



Flow Field Near a Surface 

Having found that a forcc-dipole flow describes the 
measured fiow around E. coli with good accuracy in the 
bulk (far from boundaries) , we investigated whether this 
approximation is also valid when E. coli swim close to a 
wall. Focusing 2 fim below the top of the sample cham- 
ber, and applying the same measurement technique as 
before, we obtained the fiow field shown in Fig. IE. This 
flow decays much faster than that in the bulk due to the 
proximity of a no-slip surface (Fig. IH), and the inward 
and outward streamlines are now joined to produce loops 
(Fig. IE). However, both of these differences are consis- 
tent with a simple force dipole model, and are therefore 
not due to a change in bacterial behavior. In particu- 
lar, closed streamlines are known to be a rather general 
feature of point singularities near no-slip surfaces [41[. 
Using the solution of a Stokeslet near a wall [3l| to ob- 
tain that of a force dipole near a wall yields streamlines 
(Fig. IF), and a decay (Fig. IH) of the flow field that is 
consistent with the data. The best-fit force dipole in this 
geometry yields F = 0.43 pN and £ = 2.2 ^m, consistent 
with the values obtained far from walls, but less accurate 



than those values due to the much faster decay of u{r) 
near a wall. 



Spectral Flow Analysis 

To analyze systematically the angular flow structure 
even at distances r < 6 /im, where the force dipole model 
overestimates the fiow magnitudes, it is useful to decom- 
pose the flow fleld into vector spherical harmonics. The 
resulting spectra are useful 'fingerprints' of the fiow field 
that can be compared among many different organisms, 
and against theoretical models. Such a spectral analysis 
is described in the SI Text. 



Rotational Diffusion 

Even bacteria that do not display tumbles, such as 
those studied here, do not swim in completely straight 
lines. Random changes in swimming direction due to 
thermal effects and intrinsic noise in the swimming ap- 
paratus lead to rotational diffusion which can be char- 
acterized by a coefficient Dr- From the swimming data 
recorded for E. coli far from surfaces, we measured Dr = 
0.057 rad^/s (see Materials and Methods). Even organ- 
isms that are too large to have significant thermal rota- 
tional diffusion, such as the 10-/im sized alga Chlamy- 
domonas reinhardtii, can have a significant Dr due to 
noise in the swimming mechanism |42| . From swimming 
data previously recorded for Chlam,ydom,onas [37| . we 
found Dr = 0.4 rad^/s. 



DISCUSSION 

Our measurements show that, independently of 
whether E. coli swim near or far from a surface, their fiow 
field can be described to good accuracy by a simple force 
dipole model whose parameters we determined. We now 
proceed to discuss the implications of this fiow field for 
cell-cell and cell-surface interactions. Based on the mea- 
sured parameters and the force dipole approximation, we 
calculate the effect of long-range hydrodynamics for these 
two scattering phenomena and evaluate the importance 
of fluctuations in the swimming direction. 



Hydrodynamics vs. Rotational Diffusion in Cell-Cell 
Scattering 

Fluid-mediated long-range interactions are thought to 
play an important role in collective motion in bacte- 
rial suspensions [17H23[ . These deterministic interactions 
however compete with rotational diffusion of the swim- 
ming direction. To infer the importance of longe-range 



hydrodynamics in the bulk, we consider the change in 
swimming direction of a bacterium due to hydrodynamic 
scattering with another bacterium. This can be done by 
fixing one bacterium at the origin and examining the tra- 
jectory of the other. The flow field around the bacterium 
at the origin is approximated by that of a point force 
dipole [31 1 
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where d is the unit vector in the swimming direction, and 
r is now the distance vector relative to the center of the 
dipole. The evolution of the position x and swimming 
direction d of the second swimmer in this field obeys |43j 

X = Vod + u, (2) 

d ^ -uxd + r d- E-{I ~dd), (3) 

where / is identity matrix, and the central swimmer leads 
to an advective flow it, a vorticity w, and a strain-rate 
tensor E at the position x (see SI Text for exact expres- 
sions of these quantities) . By examining the evolution of 
d{t) in a scattering process that begins at —t/2, reaches 
a minimal encounter distance r at t = 0, and ends at t/2, 
the mean squared angular change of orientation during a 
time interval t can be computed as 

(A0(t, r)^)H = (arccos[d(-i/2) • d{t/2)]^)^ , (4) 

where { ■ )h indicates an average over all possible orien- 
tations and positions of encounters. Assuming that the 
interaction time scale of the two bacteria r is small, and 
using the force dipole model from Eq. ([IJ we obtain (see 
derivation in SI) 
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where F ^ 1 is a geometric factor for E. coli. Intuitively, 
this form arises from the fact that A(/) ^ cjt, where 
the vorticity magnitude a; falls off as A/r^. To evalu- 
ate the importance of hydrodynamic interactions rela- 
tive to random fluctuations, we compare (A(/)(r, r)^)/f 
with the angular diffusion due to Brownian motion and 
intrinsic swimming variability in three dimensions, given 
by (A0(t)'^)£) = ADj-t. Balancing the effects of hydrody- 
namics and noise. 



(A0(T,rff)2)ff = (A0(r)2)c, 



(6) 



defines an effective hydrodynamic horizon r//, beyond 
which noise dominates over hydrodynamics. For scat- 
tering events with closest encounter distances r > rn 
hydrodynamics is therefore practically irrelevant. From 
this definition of rn and Eq. (O, we find 
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Due to the r^/^-dependence, the hydrodynamic horizon 
rn is rather insensitive to the particular value used for 
the interaction time scale r and, similarly, to changes 
in the other parameters. Using the measured values 
Dr = 0.057 radVs, Vq = 22 /im/s, A = 31.8 fim^/s, 
and adopting r = a/Vo — 0.1 s, where a = 3 /.jm is the 
length of the cell body, we obtain r// ~ 3.3 fim for E. coli. 
This value of th is an upper bound, as the dipolar flow 
model overestimates it for r < 6 /im (Fig. ID). At such 
small separations, however, steric repulsion, flagellar in- 
tertwining, and lubrication forces dominate the physical 
cell-cell interactions [4J, |45| . For the mean distance be- 
tween E. coli to reach rn, the volume fraction needs to 
be as high as 5 — 10%. Using our measured parameters 
in a recent theoretical calculation [461 of the critical vol- 
ume fraction for the onset of collective swimming due to 
hydrodynamic interactions leads to an even higher value, 
implying that a complete analysis of collective behavior 
in bacterial suspensions must account for steric and near- 
field interactions. 

More generally, we expect similar results to hold for 
various types of swimming bacteria (e.g., Pseudomonas 
aeruginosa, Vibrio cholerae, Salmonella typhimurium) , 
as the parameters {a,T, Dr, Fi) are similar across many 
genera. Larger organisms may display even stronger ro- 
tational diffusion due to enhanced intrinsic swimming 
stochasticity J42l |. For example, the alga Chlamydomonas 
(a ~ 10 Aim, Vo ~ 100 /im/s) has A- =i 0.4 radVs. Al- 
though the flow topology around Chlamydomonas is more 
complicated than that around bacteria, Chlamydomonas 
still produces a 1/r^ field 37[, so that our previously 



calculated result for the bacterial hydrodynamic horizon 
may be used to give an estimate oi rn ^ 7.5 /im. again 
on the scale of the organism. Thus, collisions, rather 
than long-range hydrodynamics, can also govern scatter- 
ing events of higher microorganisms. However, for organ- 
isms that produce fast fiows, have a long interaction time 
T, and a negligible Dr, like the alga Volvox, long-range 
hydrodynamic interactions are significant 47 1 . 



(7) 



Cell-Surface Scattering 

The accumulation of bacteria near surfaces is a key 
step during the initial stages of biofilm formation, and 
it has been suggested that long-range hydrodynamics 

ts an important role in bacteria-surface interactions 
As our measurements show that the field has a 
small amplitude, it is plausible that bacterial cell-surface 
scattering events could be described instead by nearly 
straight-line swimming interrupted by collisions with the 
wall which lead to alignment with the surface due to near- 
field lubrication and/or steric forces during the collision 
[28|. Our experimental results establish the key micro- 
scopic parameters required for a systematic investigation 
of whether long-range hydrodynamic interactions are rel- 



evant to bacteria-surface scattering. 

As a first step in our analysis, we performed numerical 
scattering studies by simulating the deterministic equa- 
tions of motion for an E. coli-\ike pusher force-dipole 
swimmer in the presence of an infinite no-slip surface. 
The equations of motion of the swimmer position x and 
unit orientation vector d are simply those used above 
(Eqs. (21) and ([3])) except that now u, a;, and E are 
quantities arising from the hydrodynamic image system 
in the wall (exact expressions for these quantities are 
given in the SI Text). By restricting ourselves to sim- 
ulations of the deterministic dynamics at this stage, we 
overestimate the relevance of hydrodynamic long-range 
interactions between the swimmer and the wall, as ro- 
tational diffusion of the swimming direction further di- 
minishes hydrodynamic effects. However, even without 
rotational diffusion our simulations show that long-range 
interactions of swimmers with the wall have little effect 
on the swimming dynamics (Fig. [2]). The trajectories of 
force-dipole swimmers that swim towards the wall from 
different initial angles 6o arc depicted in Fig. [5JV. The 
initial distance was chosen such that the swimmer would 
reach the wall plane {y = 0) after 1 s if hydrodynamic 
interactions were absent. Each simulation is stopped 
when a volume around the swimmer (a bacterial shape of 
length 3 ^m and diameter 0.8 ^m) crosses the wall. Fig- 
ure [23 displays the impact angle ^hit as a function of the 
initial angle 0o, illustrating that the difference between 
incidence and collision angles is small unless the swimmer 
already has a small angle of incidence. These simulations 
indicate that hydrodynamic long-range interactions are 
not likely to play an important role in cell-surface scat- 
tering for E. coll. As the swimming parameters are simi- 
lar for many bacterial genera, we again expect this result 
to apply more generally. 



Trapping by Surfaces 

When E. coli swim very close to a surface (~ 1 — 3 ^m) , 
we observed that individual bacteria spend an average 
of 64 s (standard error 4 s) next to the wall. Effective 
trapping by electrostatic attraction is unlikely, since both 
the E. coli outer cell wall and the chamber walls (bovine 
serum albumin coated onto PDMS [48[) are negatively 
charged in our liquid medium (we observed similarly long 
residence times on simple glass surfaces). However, the 
surprisingly long residence times could be caused by the 
suppression of rotational diffusion due to geometric con- 
straints on the orientation of the cell body and flagella 
near a surface. Although we showed in the previous 
section that hydrodynamics has a very small effect on 
the swimming direction before collisions with the sur- 
face, hydrodynamic attraction by the surface 127| could 
contribute to the observed trapping periods when a bac- 
terium is already very close to the surface. Considering 
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FIG. 2: Simulated dynamics of an E. co/i-like force-dipole 
swimmer near a wall. (A) Deterministic swimming trajec- 
tories towards a wall at y = 0, numerically simulated from 
Eqs. ^ and ((SI, wiiere u, a), and E are due to the hy- 
drodynamic linage system. Simulations used a time step 
At = lO^^s and the experimentally determined parameters 
A = 31.8 ^m^/s, Vo = 22 ^m/s and V = 0.88 for the force- 
dipole swimmer. The initial distance Is chosen such that the 
swimmer would reach the wall after 1 s if hydrodynamic Inter- 
actions were not present. (B) Incidence angle ^o 'vs. collision 
angle Shit with the wall for the trajectories in panel A, us- 
ing the same symbols and colors. The dotted line Indicates 
Shit = So. Both panels Illustrate that hydrodynamic long- 
range Interactions can be regarded as small perturbations for 
typical wall scattering events of E. coli. 



only hydrodynamic attraction counteracted by rotational 
diffusion, we now derive approximate expressions for the 
mean escape time t^. and escape height above the sur- 
face hf,, by mapping the underlying escape process onto 
a Kramers' problem [49|,|50| for the noise-induced escape 
over a potential barrier. The main arguments and impli- 
cations are summarized below, while a detailed derivation 
is given in the SI Text. 

We again approximate the E. coli flow field by the 
dipole model, because a force-dipole placed close to a wall 
accurately captures the measured flow field parallel to the 
surface (see Fig. lE-H). Thus, EclJTT]) is modified to ac- 
count for the presence of the wall [3l| , as discussed in the 
SI Text, and near-field hydrodynamic lubrication effects 
are neglected. A bacterium is able to escape from the 
surface, if its swimming velocity component perpendic- 
ular to the surface, Vosin0, exceeds the attraction from 
its hydrodynamic image (see SI Text), which yields the 
defining relation for the escape angle 9e, 



sin6'e = A[l-3(sin6'e)^] , A 



3A 
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(8) 



For E. coli swimming at distances h > 1.5 /im from the 
wall, the escape angles are small, Oe{h) < 11° ^ 1 rad 
so that linearization of Eq. ([S]) is a good approximation, 
giving 9e ^ A. 

After colliding with the wall, a bacterium may have a 
small positive angle 9 < de with the surface. The equa- 
tion of motion for 9 can then be rewritten as a Langcvin 
equation [49l | in terms of the derivative of an effective 'po- 



tential' U{9), and a diffusion term with Gaussian white 
noise ^(t), 



-^ + (2z?:)V2e(i), 



^(^) - 2;: 



(9) 



where the approximation 9 <^ 1 reduces U{9) to a har- 
monic potential, yielding a time scale k = 16/i'^/(9A) 
that characterizes hydrodynamic realignment. D* is the 
rotational diffusion constant close to the surface in the 
direction perpendicular to the surface, which is expected 
to be smaller than our measured value Dr = 0.057 rad^/s 
far from boundaries, due to geometric constraints on the 
bacterial orientation near a surface. The generic form 
of this Langevin equation means that finding the res- 
idence time for a bacterium near a wall is a Kramers 
problem [49|, |50| for the escape over a barrier AC/. As 
the organism can escape if 6* > 9e, we have AU = U{9e)- 
By considering the height at which AC/ = D*, i.e., the 
distance at which hydrodynamic torque is comparable 
to diffusion, we can obtain an expression for the escape 
height 



_ 1 /81 A^ 
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(10) 



Using our measured values for E. coli, we find he = 
1.7 /xm X {Dr/D*)^/"^ , illustrating that hydrodynamics is 
practially negligible if E. coli are more than a cell length 
away from the wall. 

As long as the torque exerted by the hydrodynamic im- 
age is small, AC/ ^ D*, the typical escape time is set by 
the rotational diffusion time scale 0^/D*. For high bar- 
riers AC/ » D* (in practice, AC/ > ?>D* often suffices, 
yiel ding /i < 1.5 x {Dr/ 0*^/"^ IJ-i^), transition state the- 
ory [49[ implies that the mean escape time is modified by 
an Arrhenius-Kramers factor, so that approximately 



te{h) 
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(11) 



Using the quadratic approximation for AC/, and Eq. 
(ITU|) to express D* in terms of he, we find that tg oc 
cxp{he/hy. This dramatic scaling arises from the fact 
that the dipole model overestimates the flow field close to 
the bacterium, but generally hints at the possibility of a 
strong, hydrodynamically induced increase of ie when the 
cells get closer to the surface. We may also evaluate Eq. 
([TT|) at a height h = 1.5 /xm, where both the Arrhenius- 
Kramers factor and the approximation 9e{h) <^ 1 are 
valid, to give 



t, w0.78sx 



exp 
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The latter estimate suggests that hydrodynamic effects 
can possibly explain the experimentally observed long 
residence times near a wall, even for values of Dr/ D* 



that are only moderately larger than 1. It is, however, 
important to note that this expression for t^ presents an 
upper bound, because the dipole model overestimates the 
actual flow flcld at distances < 6 /im from the bacterium 
(Fig. IH), even though the model still correctly captures 
the flow topology. 

The above considerations show that hydrodynamics 
is negligible if a bacterium is more than a body length 
away from the wall, but that hydrodynamic effects may 
contribute to the experimentally observed long residence 
times of bacteria close to no-slip surfaces. A more de- 
tailed understanding of the escape problem remains an 
important future challenge, requiring new methods for 
measuring D* and further theoretical studies of the near- 
field interactions between bacteria, their flagella, and sur- 
faces. However, even if a more accurate description of the 
hydrodynamics should become available in the future, 
one can still expect the mean escape time to follow an 
Arrhenius-Kramers law of the form (|lip with a suitably 
adapted effective potential U and additional prefactors 
that account for the curvature at the barrier 14911. 



CONCLUSIONS 

We have presented the first direct measurement of 
the flow field generated by individual freely- swimming 
bacteria, both in bulk fluid and close to a solid sur- 
face. For distances > 6 /im, the experimentally measured 
flow field is well-approximated by a force-dipole model; 
at smaller distances the dipole model overestimates the 
flow. Generally, the flow field of E. coli differs markedly 
from those created by higher microorganisms, such as 
Chlamydomonas |37l.l38| and Volvox 37 1. With regard to 
the future classification of flow flelds of microorganisms, 
a decomposition in terms of vector spherical harmonics 
can provide a useful systematic framework, similar to the 
classiflcation of the electronic orbital structures in atoms 
or molecules. 

Theories of collective behavior in bacterial suspensions 
identify as a fundamental process the pairwise interaction 
of bacteria, often assumed to be dominated by long-range 
fluid flows established by the action of swimming [25|. 
Our analysis suggests that noise, due to orientational 
Brownian motion and intrinsic swimming stochasticity, 
drowns out hydrodynamic effects between two bacteria 
beyond a surprisingly small length scale of a few mi- 
crons. This implies that hydrodynamic effects will be 
relevant only in sufficiently dense bacterial suspensions. 
However, under such conditions, the flow structure close 
to the bacterial body and contact interactions (e.g., flag- 
ellar bundling, stcric repulsion) will be more important 
than the asymptotic long-range details of individual mi- 
croswimmer flow fields. 

Insights into the biochemical and physical interactions 
between bacteria and surfaces are crucial for understand- 



ing the dynamics of biofilm formation, the emergence 
of collective bacterial behavior in boundary layers and, 
thus, more generally the evolution from unicellular to 
multicellular, cooperative forms of life. Our results sug- 
gest that long-range hydrodynamic effects play a negligi- 
ble role in the scattering of E. coli with surfaces before 
collisions. However, hydrodynamic effects can, at least 
partially, account for the observed trapping of bacteria 
within a few microns of the surface. The analysis pre- 
sented herein lends support to the hypothesis |5l| that 
turbulent swarming patterns in bacterial films arise pri- 
marily due to steric repulsion and other near field inter- 
actions, since long-range hydrodynamic interactions be- 
come suppressed near surfaces due to cancellation effects 
from the image swimmer. 

Our experimental and theoretical results favor 
collision-dominated models |28l43Cl| for the accumulation 
of bacteria at surfaces over models based on long-range 
hydrodynamics J27l | . To obtain a more complete dynam- 
ical picture of biofilm formation, future efforts should 
focus on developing more precise measurement meth- 
ods and advanced models that include lubrication effects 
and biochemical bacteria-surface interactions. While our 
combination of measurements, simulations, and theory 
shows that long-range physical interactions are negligible 
for bacterial cell-surface scattering, fluid-mediated cou- 
pling could become important for organisms swimming 
against or in contact with a surface, as the organism 
is then no longer force-free, resulting in a substantially 
longer range of hydrodynamic interactions [47|, |52| . 

However, the main implication of the present study is 
that short-range forces are likely to dominate the interac- 
tions between swimming bacteria, so that collective mo- 
tion in bacterial suspensions, thin films [J, |53[ and thin 
wetting layers [SJ] relates closely to that seen in driven 
granular systems ISSl. a ssemblages of biofilaments [56[, 



and animal flocks (57|,|58[. This suggests that many of the 
principles that determine flocking and self-organization 
in higher animals should also govern the collective mo- 
tion of the smallest organisms. 



broth (1% tryptonc, 0.5% NaCl) and 100 A^g/mL ampi- 
cillin. A single-colony isolate from an overnight plate was 
used to innoculate 10 mL of T-broth containing ampi- 
cillin and 0.1 mM isopropyl /3-D-l-thiogalactopyranoside 
(IPTG, Sigma), which was then grown for 7 h on a rotary 
shaker (200 rpm) at 33 °C. This culture was diluted 1:1 
with fresh T-broth containing ampicillin and IPTG as 
above, 0.2% bovine serum albumin, and 0.2 ^m fluores- 
cent microspheres (505/515, F8811, Invitrogen) at con- 
centration 9 X 10^ beads/mL. This bacterial suspension 
(~ 1.6 X lO'' cells/mL) was loaded into a polydimethyl- 
siloxane (PDMS) microfluidic device consisting of cylin- 
drical measurement chambers (height 100 /im, radius 750 
/im) connected by thin channels. After filling the device, 
it was sealed to reduce background fluid motion. 



Measurement of the Flow Field 

Using a Zeiss Axiovert inverted microscope with a 40 x 
oil objective (NA 1.3), we simultaneously imaged bacte- 
ria and microspheres under fluorescence conditions at 40 
fps (Pike, Allied Vision Technologies) and at a temper- 
ature of 24 ± 1 °C. To measure the flow field far from 
walls, we focused on a plane 50 /im inside the chamber 
to minimize surface effects. To measure the flow field 
close to a no-slip surface, we focused on a plane 2 fim be- 
low the top surface of the sample chamber. Each movie 
was analyzed with custom Matlab software that precisely 
tracked bacteria by fitting an ellipsoidal two-dimensional 
Gaussian shape. For each cell swimming along the focal 
plane for > 1.5 s, we collected the instantaneous veloc- 
ity of all fluorescent tracers up to a distance of 75 /im, 
using standard particle tracking algorithms. The result- 
ing '--^ 5 X 10^ tracer velocity vectors were binned into a 
0.63 /im square grid (shown in Fig. lA, E). The mean of 
the well-resolved Gaussian distribution in each bin was 
taken as a local measure of the flow field. To measure 
the mean residence time of bacteria near a surface, we 
used the movies that were recorded for measuring the 
flow field near the wall. 



MATERIALS AND METHODS 

A detailed description of the mathematical models is 
provided in the SI Text. The experiments are summa- 
rized below. 



Culture Conditions 

We used Escherichia coli strain HCB437 carrying the 
plasmid pEGFP (Clontech, BD Biosciences) , kindly sup- 
plied by Douglas B. Weibel (University of Wisconsin- 
Madison) and Howard C. Berg (Harvard University). 
Cells were streaked on 1.5% agar plates containing T- 



Measurement of the Rotational Diffusion 

From the tracks of E. coli that swam in the focal plane 
for > 1.5 s, at a distance of 50 /im from the top and 
bottom surfaces, we determined an average swimming 
direction at time t by using the direction between the 
bacterial positions at f — 0.05 s and t + 0.05 s. Com- 
puting the change in average swimming direction A(j) 
revealed diffusive scaling, so that we obtained Dr from 
the equation for two-dimensional orientational diffusion, 
(|A(/)p) = 2DrAt, over a time interval At. We measured 
Dr for Chlamydomonas with the same procedure, using 
cell-tracking data described earlier [37[ . 
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